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Abstract 

The problem of the stability of a nonlinear thermomagnetic wave with respect to small thermal 
and electromagnetic perturbations in hard superconductors was studied. It is shown that spatially 
bounded solutions may correspond only to the perturbations decaying with time, which implies 
stability of the nonlinear thermomagnetic wave. 

Previously [1], it was demonstrated that, depending on the surface conditions, stationary 
nonlinear thermomagnetic waves of two types, E and H, may exist in a superconductor. 
However, the problem of stability of the nonlinear waves with respect to small thermal and 
electromagnetic perturbations in superconductors is insufficiently studied. We will consider 
the stability of a nonlinear thermomagnetic wave with respect to small thermal and elec- 
tromagnetic perturbations in a superconductor occurring in the critical state. It will be 
demonstrated that spatially bounded solutions may correspond only to the perturbations 
decaying with time. The evolution of the thermal (T) and electromagnetic (E and H) per- 
turbations in a superconductor is described by a nonlinear thermal conductivity equation 
[2] 

u(T)^ = A[ K (T)AT]+jE, (1) 



in combination with the Maxwell equations 

l dH 

c dt ' 



-* ldH -> 4tt^ 

rotE = — , rotH = — j (2) 



and the critical state equation 

j = j c (T,H)+j r (E). (3) 

where v = u(T) and k = k(T) are the heat capacity and thermal conductivity coefficients, 
respectively; j c is the critical current density; and j r is the resistive current density. Let 
us consider a plane semiinfinite sample (x > 0), exposed to an external magnetic field 
-» dH 

H = (0,0, H e ), increasing with the time at a constant rate = 0. According to the 

dt 

Maxwell equation (2), the sample will feature a vortex electric field E = (0,.E e ,0) parallel 
to the current density vector: j: E \\ j. 

In order to describe the j c {T, H) function, we will use the Bean-London equation of the 

dj c 

critical state (-^ = 0) [3]. According to this model, the j c (T) function can be presented as 
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j(T) = Jc (T ) 



T-T 



T c — T n 



: where j = j c (T ) is the equilibrium current density, T c is the 

critical temperature, and To is the cooler temperature. The characteristic form of j r {E) in 
the region of sufficiently strong electric fields (E > Ef) can be approximated by a piecewise 
linear function j r m (JfE. In the region of small field strengths E < Ef the j r (E) function 
exhibits a nonlinear character related to a thermoactivated flux creep [4] . For an automodel 
solution of the type £ = x — v t describing a wave propagating at a constant velocity v along 
the x axis, the system of Eqs.( 1 )-(3) acquires the following form: 

-v[N(T)-N(T )] = ^-^-E 2 , (4) 



dE Anv . 



E = -H. (6) 

c 

where N(T) = ]v{T)dT. 
o 

Excluding variables T and H with the aid of expressions (4) and (6) and taking into 
account the boundary condition E(z — > — oo) = E e , we arrive at a differential equation 
describing the distribution E(z) 

d 2 E dE 4ttv 2 E 2 

Here z — -7- and (3 = are dimensionless parameters, L = - — - the magnetic field 
L L A-kjq 

penetration depth into the superconductor, r = — — is the ratio of the coefficients of thermal 

Dm 

Kj C 2 U Ij 2 

D t — — and magnetic D m = diffusion, t K = is the thermal diffusion time, and 

V A-KGf K 

K 

E K = — - is a constant parameter. 
aL 2 

In the approximation of a weakly heated superconductor (T — T ) << T , the heat capacity 
v and thermal conductivity n coefficients are weakly dependent on the temperature profile. 
As is well known (see, e.g., [4]), the magnetic flux variations in a hard superconductor 
occur at a much greater rate as compared to those of the heat transfer, so that r << 1 or 
Dt « D m . In this approximation, we may neglect the terms related to dissipative effects 
in Eq. (7), after which a solution to this equation can be presented in the following form [5] 



(3t 

1 - th—(z - zq) 



(8) 



Relationship (8) describes the profile of a thermomagnetic shock wave propagating into the 
superconductor. In order to study the stability of a nonlinear wave with respect to small 
perturbations, it is convenient to write a solution to Eqs. (l)-(3) in the following form: 
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T(z,t) = T(z) + ST(z,t)exp 



Xt 



E(z, t) = E(z) + 5E(z, t) exp 



(9) 



where T(z) and E(z) are the stationary solutions and ST, SE are small perturbations. 
Substituting expressions (9) into Eqs. (1)- (3), assuming ST, 5E « T,E in the limit of 

t « 1 (this corresponds to a "fast" instability, A >> [4]), we obtain an equation for 

K 

determining the eigenvalues of A [6] 

Using the substitution of variables of the type 

— 

£ = thy, 1 - £ = 2s, and e(£) = (1 - £ ) 2 $(0, Eq. (10) can be reduced to a standart 
hyp ergeo metric equation 

< l ~ ~ [ Q + 1 - 2s ( fi + !)]^ " " !) " P(P " = ( U ) 

Here, the even and odd integrals can be presented as 

$i = F(e-l,e + 2,e+l,s), (12) 

$ 2 = s - e F(-l,2,l-e,s), (13) 

where F is the hypergeometric function [7]. The function must be finite at the singular 
point s — 1. The e values for which $x is finite evidently correspond to a discrete spectrum: 
6i — 0, 1 or Aj = — r _1 ,0. The function $ 2 is finite only for e = 0. Any solution in the 
form of a running wave is char- acterized by translational symmetry, which implies that a 
"perturbed" stationary profile E(z) determined by Eq. (8) corresponds to the eigenvalue 

of the ground state A = 0. Differentiating Eq. (7) with respect to z, one may readily see 

dE dE 
that — - — is an eigenfunction corresponding to A = 0. Indeed, the perturbation 5E = — — 

dz dz 
essentially represents a small wave displacement. Thus, we may suggest that tha function 

dE 

— — exponentially tends to zero for z — > +oo and the corresponding eigenvalue is zero. 
Therefore, the problem cannot possess positive eigenvalues and ReXi < 0. This result implies 
that the wave is stable with respect to relatively small thermal ST and electromagnetic SE 
fluctuations. The analysis of the second linearly-independent solution leads to the same 
conclusion. 

REFERENCE 



[1] I.L.Maksirnov,Yu.N.Mastakov, and N.A.Taylanov, Fiz. Tverd. Tela (Leningrad) 28 (8), 2323 (1986) [Sov. Phys. 
Solid State 28, 1300 (1986)]. 



3 



[2] N.A.Taylanov Superconduct. Science and Technology, 14, 326 (2001). 

[3] C. P. Bean, Phys. Rev. Lett. 8 (6) ,250 (1962). 

[4] R. G. Mints andA. L. Rakhrnanov, Instabilities in Superconductors (Nauka, Moscow, 1984). 

[5] . V. I. Karprnan, Nonlinear Waves in Dispersive Media (Nauka, Moscow, 1973; Pergarnon, Oxford, 1975). 

[6] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Nonrelativistic Theory (Fizrnatgiz, Moscow, 1963; Per- 
garnon, New York, 1977). 

[7] D. S. Kuznetsov, Special Functions (Vysshaya Shkola, Moscow, 1965). 



4 



